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The exponential sensitivity of cluster number counts to the properties of the dark energy implies 
a comparable sensitivity to not only the mean but also the actual distribution of an observable mass 
proxy given the true cluster mass. For example a 25% scatter in mass can provide a ~ 50% change 
in the number counts at z ~ 2 for the upcoming SPT survey. Uncertainty in the scatter of this 
amount would degrade dark energy constraints to uninteresting levels. Given the shape of the actual 
mass function, the properties of the distribution may be internally monitored by the shape of the 
observable mass function. An arbitrary evolution of the scatter of a mass-independent Gaussian 
distribution may be self-calibrated to allow a measurement of the dark energy equation of state of 
cr(w) ~ 0.1. External constraints on the mass variance of the distribution that are more accurate 
than < 0.01 at 2 ; ^ 1 can further improve constraints by up to a factor of 2. More generally, 

cluster counts and their sample variance measured as a function of the observable provide internal 
consistency checks on the assumed form of the observable-mass distribution that will protect against 
misinterpretation of the dark energy constraints. 


I. INTRODUCTION 

It is well known that cluster counts as a function of 
their mass are exponentially sensitive to the amplitude 
of the linear density field and hence the dark energy de¬ 
pendent growth of structure. Unfortunately the mass of 
a cluster is not a direct observable and their numbers can 
only be counted as a function of some observable proxy 
for mass. Typical proxies include the Sunyaev-Zel’dovich 
flux decrement, X-ray temperature, X-ray surface bright¬ 
ness or gas mass, optical galaxy richness, and the weak 
lensing shear. The exponential sensitivity to mass trans¬ 
lates into a comparable sensitivity to the whole distribu¬ 
tion of the observable given the mass not just the mean 
relationship. 

While scatter in the observable-mass relation is typi¬ 
cally addressed in studies of the local cluster abundance 
(e.g. Q), it is commonly ignored in forecasts for upcom¬ 
ing high redshift surveys (e.g. Hill) While it is true 
that scatter in the observable of a known form does lit¬ 
tle to degrade the dark energy information, uncertainties 
in the distribution directly translate into uncertainties in 
the dark energy inferences that must be controlled. 

In this Paper we undertake a general study of the im¬ 
pact of uncertainty in the observable-mass distribution 
on high redshift cluster counts. Previous work on fore¬ 
casting prospects for dark energy constraints have exam¬ 
ined the effect of scatter under specific and typically more 
restrictive assumptions. For example, the change in the 
number counts, known as Eddington bias [3, has been 
assessed for a fixed cut in signal-to-noise of cluster detec¬ 
tion via the Sunyaev-Zel’dovich flux in a hydrodynamic 
simulation and through modeling a constant scatter 
in mass 0. However it is the uncertainty in the scatter, 


or the error in the correction of the bias, that degrades 
dark energy constraints. Along these lines, Levine et al. 
Q considered the marginalization of a constant scatter 
in the mass-temperature relation for clusters but with 
strong external priors on the dark energy parameters. 

Prospects for the self-calibration of the mean 
observable-mass relation have been extensively studied 
recently. Self-calibration relies on the fact that both the 
shape of the mass function and the clustering of clus¬ 
ters [T^ can be predicted from cosmological simulations. 
Much of the information in the latter can be extracted 
from the angular variance of the counts so that costly 
spectroscopy can be avoided Thus by demanding 

consistency between the counts and their sample vari¬ 
ance across the sky as a function of the observable mass, 
one can jointly solve for the cosmology and the mean ob¬ 
servable mass relation. Here we show that the shape of 
the mass function is even more effective at monitoring 
the scatter in the observable-mass relation. 

We begin in JUlwith a discussion of our parameteriza¬ 
tion of the observable-mass distribution and assess its im¬ 
plications for the dark energy. We examine the prospects 
for self-calibration in m and conclude in 


II. OBSERVABLE-MASS DISTRIBUTION 

The cosmological utility of cluster number counts 
arises from their exponential sensitivity to the amplitude 
of the linear density field. For illustrative purposes, we 
will employ a fit to simulations for the mass function or 
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FIG. 1: Scatter of crinM = 0.25 in the observable-mass rela¬ 
tion changes the mass distribution of clusters above an observable 
threshold (curves) to provide an excess of clusters scattering 

up (dark shaded) versus down (light shaded) across the threshold. 
Here the intrinsic mass function (thick line) has been normalized 
to Mq = Mq and evaluated at 2 : = 0. 


the differential comoving density of clusters 




dn 

d\nM 




FIG. 2: Local power law index a of the mass function as a function 
of mass dn/dlnM oc for 2 = 0, 0.5 and 1. The relative impor¬ 
tance of scatter versus bias can be scaled through a and Eqn. e 
to alternate mass and redshift ranges than considered here. 
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- [erfc(xi) - erfc(a::i+i)] , (4) 


where a^{M] z) = aj^{z), the linear density field variance 
in a region enclosing M = A-kR^ pm/i at the mean matter 
density today pm- 

To exploit this exponential sensitivity, the observable- 
mass distribution must be known to a comparable accu¬ 
racy. Let us consider the probability of assigning a mass 
M°'^® to a cluster of true mass M to be given by a Gaus¬ 
sian distribution in InM as motivated by the observed 
scatter in the scaling relations between typical cluster 
observables (e.g. Q) 

p{M°'^^\M) = exp , (2) 

V 


where 




In M°’^" - In M - In 

\/2^1n M 


( 3 ) 


For simplicity we will allow the mass variance ^^.d 

the mass bias In to vary with redshift but not mass. 

We implicitly exclude sources of scatter due to noise in 
the measurement of M°'°^ which certainly would depend 
on M°'°^ but in a way that is known given the proper¬ 
ties of a specific survey. More generally, our qualitative 
results will hold so long as any trend in mass at a fixed 
redshift is known. 

The average number density of clusters within a range 
defined by cuts in the observable mass Mf°^ < M°'°^ < 


where Xi = x{Mf°^). The mean number of clusters in a 
given volume Vi is then 

fhi = niVi. (5) 

Note that in the limit that > 0 and M°'^\ —> cx), 

rhi is the usual cumulative number counts above some 
sharp mass threshold. 

An unknown scatter or more generally uncertainty in 
the distribution of the observable mass given the true 
mass causes ambiguities in the interpretation of num¬ 
ber counts. Fig. n shows the expected mass distribu¬ 
tion of clusters above a certain M°'°^ given a scatter of 
trinM = 0.25. As the observable threshold reaches the 
exponential tail of the intrinsic distribution, the excess 
of upscattered versus downscattered clusters can become 
a significant fraction of the total. Since at high redshift 
a fixed M°'^® will be further on the exponential tail, even 
a constant but unknown scatter can introduce a trend in 
redshift that will degrade the dark energy information in 
the counts. 

The relative importance of scatter can be understood 
by examining the sensitivity of the counts to the scatter 
around m ~ ^ 

dhifhi 1 

hm 5-= — —-5— 

2n,dlnM 

Thus the steepness of the mass function around the 
thresholds in the observable mass determines the excess 
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FIG. 3: Fractional sensitivity of the number counts in redshift 
to the mass variance scaled to (0.25)^ and 

0.25 respectively. The relative importance of the variance increases 
with 2 . 


due to upscatters. Note that it is the variance 
rather than the rms scatter crinM that controls the up- 
scattering effect. For example, since 


lim 


d\nmi 9 In m 

= 2(TlnM 


dcJln 


M 


d^fnM 


= 0 , 


( 7 ) 


the sensitivity to the rms scatter depends on the true 
value of the scatter and vanishes at (Ti„ m = 0. Conversely 
from Eqn. an observable with say half the scatter 
would have a quarter of the fractional effect on number 
counts in this limit. 

On the other hand the sensitivity to the bias is given 
by 


lim 


9 In iTii 
0 d In 


1 dn Xi—O 
fiidhiM xi+i=o 


( 8 ) 


Thus the relative importance of scatter compared with 
bias can be estimated through the local power law slope 
of the mass function dn/dlnM oc (see Fig.® 


Id'^n/dln^M 1 

-= — cx(M) 

2 dfi/dhiM 2 


(9) 


Uncertainties in scatter can dominate those of bias for 
the steep mass function at high mass or redshift. 

These expectations are borne out at finite scatter by 
a direct computation of the number count sensitivity. 
Fig. m shows the sensitivity of number counts above 
Mohs ^ in redshift bins of Az = 0.1 eval¬ 

uated around InM'^^®'® = 0 and a finite scatter cti^m = 
(0.25)^. Both in terms of absolute sensitivity and rela¬ 
tive sensitivity compared to the bias, the importance of 
scatter increases with redshift. Uncertainties in the mass 
variance of = (0.25)^ would produce a ~ 50% 

uncertainty in the number counts at z = 2. For the high 



FIG. 4: Simulations predict the sample variance of counts, shown 
for a typical volume {R = 50h~^ Mpc, upper curves), and the 
shape of the mass function (lower curves) as a function of mass. 
Self-calibration is assisted by binning the selection, multiplied here 
by 0.01 for clarity and shown here with = 0.25, into 5 bins 

of Alog]^QM°^® = 0.2 (solid lines) as opposed to a single thresh¬ 
old binning of M°'^® > Mq (dotted line). With only 

threshold binning, joint changes to the cosmology, mass bias, and 
scatter are degenerate with the dark energy equation of state w 
(long dashed lines). 


z counts to provide cosmological information, the scatter 
must be known to significantly better than this level. 

Given that the relative effect of scatter depends on the 
local slope of the mass function, measuring the counts 
as a function of monitors the scatter in the mass- 

observable relation. Combined with additional informa¬ 
tion in the sample variance of the number counts, an 
unknown evolution in and ^ may be internally 

calibrated. 


III. SELF-CALIBRATION 


To assess the impact of uncertainties in the observable- 
mass distribution, we employ the usual Fisher matrix 
technique. For illustrative purposes we take a fiducial 
cluster survey with specifications similar to the planned 
South Pole Telescope (SPT) Survey: an area of 4000 deg^ 
and a sensitivity corresponding to a constant = 

We further divide the number counts into 
bins of redshift Az = 0.1 from an assumed optical photo¬ 
metric followup out to z = 2 and 400 angular cells of 10 
deg^ for assessment of the sample variance of the counts 
(see for an exploration of these choices). Finally to 
study the efficacy of self-calibration from binning of the 
observable, we compare 5 bins of Alogj^g = 0.2 ver¬ 
sus a single bin of (see Fig.®. 

The Fisher matrix is constructed out of predictions 
for the number counts and their covariance. The mean 
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FIG. 5: Efficacy of self-calibration through (a) the clustering information in the sample variance alone, (b) the shape of the counts in 5 
mass bins and (c) both (note the xlO change in scale). From outer to inner ellipse, each at the 68% CL, the assumptions on the redshift 
evolution of the bias and scatter are tightened from a free functional form of 20 parameters each in bins of /S.z = 0.1, a power law evolution 
in an additional cubic form for the evolution in Without any form of self-calibration no relevant dark energy constraints 

are possible under any of these assumptions. For reference we show the baseline results of a completely fixed and as the solid 

innermost ellipse 


number counts rrii possess a sample covariance 


of a 


Sij = ((m* - mi){mj - nij)) 


bimibjTTij 


(Pk 


Wi*(k)VF,(k)P(fc), (10) 


given a linear power spectrum P{k) and the Fourier 
transform of the selection window Wi(x). The pixel in¬ 
dex i here runs over unique cells in redshift, angle, and 
observable mass. Here bi is the average bias of the clus¬ 
ters predicted from the distribution in Eqn. 0 


dM dm 


where we take a fit to simulations of M 
acSljP-l 2pc 


b{M-z) = 1 + 


i5c[l -I- {a5l/PYc\ 


( 11 ) 


( 12 ) 


with Oc = 0.75, Pc = 0.3, and 5c = 1.69. That the sample 
variance, or the clustering of clusters, is a known function 
of mass provides a second means of self-calibration M- 
Note that for a given volume defined by the redshift and 
solid angle, the mean counts for different ranges in the 
observable mass are fully correlated. Finally the 

total covariance matrix is the sample covariance plus shot 
variance 


Cij — Sij TTliSij . ( 10 ) 

The Fisher matrix quantifies the information in the 
counts on a set of parameters Pa as [Hill 

Fcp = + iTr[C-iS,„C”iS,;3], (14) 

where the first piece represents the information from 
the mean counts and the second piece the informa¬ 
tion from the sample covariance of the counts. We 


have here arranged the counts per pixel i into a vec¬ 
tor m = (mi,... jmiVp,^) and correspondingly their co- 
variance into a matrix. The Fisher matrix approximates 
the covariance matrix of the parameters Ca /3 « [F~^]q ;3 
such that the marginalized error on a single parameter 
is Cr(Pa) = [F~'^]aa- When considering prior informa¬ 
tion on parameters of a given a{pa) we add to the Fisher 
matrix a contribution of <J~'^{Pa)5ai3 before inversion. 

For the parameters of the Fisher matrix we begin with 
six cosmological parameters: the normalization of the ini¬ 
tial curvature spectrum 5c^{= 5.07 x 10“®) at fc = 0.05 
Mpc“^ (see [il for its relationship to the more tradi¬ 
tional tTg normalization), its tilt n(= 1), the baryon den¬ 
sity Q.bh?{= 0.024), the dark matter density Hmh^(= 
0.14), and the two dark energy parameters of interest: 
its density Hde(= 0.73) relative to critical and equation 
of state w{= —1) which we assume to be constant. Val¬ 
ues in the fiducial cosmology are given in parentheses. 
The first 4 parameters have already been determined at 
the few to 10% level through the CMB and we will 
extrapolate these constraints into the future with priors 
of (T(ln(5(;) = (j{n) = CT(lnf2bh^) = aihi^mP) — 0.01. 

For the observable-mass distribution we choose a fidu¬ 
cial model of = 0 and — (0.25)^. 

The results that follow do not depend on the specific 
choice as we have explicitly shown by testing a much 
smaller fiducial scatter of — (0.05)^. 

Given an observable-mass distribution fixed at the 
fiducial model and the priors on the other cosmologi¬ 
cal parameters, the baseline errors on the dark energy 
parameters are a{VLi^,w) = (0.008,0.03). The mere pres¬ 
ence of scatter in an observable does not necessarily de¬ 
grade the cosmological information; in fact for reasonable 
scatter it actually enhances the information by effectively 
lowering the mass threshold at high redshift. 

However, cosmological parameter errors are degraded 
once observable-mass parameters are added in a joint 
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FIG. 6: Prior knowledge of the scatter in the observable-mass 
relation at the level of < (0.1)^ for each of the 20 redshift 

bins can improve dark energy constraints beyond self-calibration 
by a factor of 2 (solid line). This cumulative effect of independent 
priors is compared with the joint effect of priors on the 4 parameters 
of a a{Ba) = m') ^ cubic which captures most of 

the residual uncertainty in Fig.[3 The latter reflects uncertainties 


fit. As the most general case, we take independent 
and crinM('^i) parameters for each redshift 
bin for a total of 40 parameters. As discussed in s]nl for 
the Fisher results to be valid around a fiducial cinM Ij 
the mass variance and not its scatter must be chosen as 
the parameters. Because the evolution in the cluster pa¬ 
rameters is expected to be smooth in redshift we alterna¬ 
tively take a more restrictive power law evolution in the 
bias 


In= Ab-I-nt, ln(l-I-Zi) (15) 

and/or a Taylor expansion of around z = 0 


a 


2 

InM 


{Zi) = <7 


2 

InM 


fid 


N„-l 

+ E 


a=0 


(16) 


With no self-calibration, i.e. no clustering information 
from the sample variance and no binning in in¬ 

teresting constraints on the dark energy are not possible 
even for the restricted evolutionary forms of Eqns. m- 
(tT?^ if crinM('^) is allowed to evolve (V^ > 2). Even 
restricting the parameters to a single constant scat¬ 
ter {Na = 1) causes a degradation to a{nA,w) = 
(0.37,0.24). 

As shown in Fig. [padding in the sample (co)variance 
information in the Fisher matrix of Eqn. 03 for a sin¬ 
gle bin in helps but still does not allow for full 

self-calibration of an arbitrary evolution in (7‘f^j^{z) even 
when is restricted to power law evolution. Fur¬ 

ther restricting the evolution in the mass variance to a 
cubic form = 4 yields a{nA,w) = (0.22,0.17) and to 
a constant form = 1 yields a{il.A,'w) = (0.15,0.07). 


Employing the information contained in the shape of 
the counts through binning allows for a more robust 
self-calibration. In the case of arbitrary evolution for both 
the bias and the scatter a{flA,w) = (0.03,0.21). With 
a power law form for the bias a{flA,w) = (0.02,0.11); 
with an additional cubic form for the mass variance 
a{nA,w) = (0.02,0.10); with a constant form for the 
scatter cr(r2A,u>) = (0.02,0.06). 

External priors on the observable-mass distribution 
from simulations and cross-calibration of observables can 
further improve on self-calibration. Cross-calibration of 
cluster observables may involve a subsample of clusters 
which have detailed mass modeling from lensing or X- 
ray temperature and surface brightness profiles assum¬ 
ing hydrostatic equilibrium M- In Fig.IHlwe explore the 
effect of independent priors on the 20 parame¬ 
ters in the power law context. Priors of the level 

of cr(cr[^^) = (0.1)^ would suffice to improve a{w) by a 
factor of 2. Note that the potential further improvement 
in w errors comes from the ability to change the scatter 
smoothly from z = 0 to z ^ 1. Since we take the priors to 
be independent, their cumulative effect implicitly poses 
a much more stringent constraint on the possible smooth 
evolution of than any one individual prior. 

To better quantify the implications of the joint prior, 
note that the self-calibration errors on the cubic form 
in Fig. HI: nearly coincide with the fully arbitrary form. 
Taking independent priors on the 4 Ba parameters of 
a{Ba) = 2a{af^n^) to reflect the assumed uncertainty 
around z = 1 yields the dashed curve in Fig. HI The 
full improvement requires priors at the (0.02)^ level 
and a minimum of (0.1)^ for substantial improvements. 
If these priors are to come from mass modeling of observ¬ 
ables then a fair sample of more than ~ 100 clusters at 
z ~ 1 with accurate masses will be required. Accurate 
masses will be difficult to obtain at the low threshold of 
employed here. 

Priors on can also improve constraints. For a 

completely fixed errors for an arbitrary evolution 

in crfjjAf <7{V,a,w) = (0.01,0.06). Conversely for a 
completely fixed errors for an arbitrary evolution 

of jy^bias cr{nA,w) = (0.02,0.13). 

Finally to assess the possible impact of unknown trends 
in the mass bias and variance with mass at fixed redshift, 
we limit the bins to 2 separated by A log^Q = 
0.2 from the threshold. In this case errors for arbitrary 
evolution degrade slightly to cr(r2A,w) = (0.03,0.26) 
and with power law mass bias and cubic variance to 
a{V,A,w) = (0.03,0.11). Thus most of the informa¬ 
tion from self-calibration comes from the small range in 
masses around the threshold reflecting the steepness of 
the mass function. The mass bias and mass variance 
need only be constant or slowly varying in a known way 
in mass across a range in masses comparable to the ex¬ 
pected scatter for self-calibration to be effective. In any 
case, bins at higher masses also monitor the validity of 
this assumption in practice. 
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IV. DISCUSSION 

The exponential sensitivity of number counts to 
the cluster mass requires a calibration of the whole 
observable-mass distribution before cosmological infor¬ 
mation on the dark energy can be extracted. We have 
shown that even in the case of an unknown arbitrary 
evolution in the mass bias and scatter of a Gaussian dis¬ 
tribution there is enough information in the ratio of the 
counts in bins of the observable mass and their sample 
variance to calibrate the distribution and provide inter¬ 
esting constraints on the dark energy. 

For the more realistic case of an unknown power law 
evolution in the mass bias the forecasted errors for 

the fiducial SPT-like survey are a{fl/y,w) = (0.02,0.11). 
To further improve on these constraints, external con¬ 
straints on the mass variance would need to achieve an 
accuracy of cr(cr[^^) < (0.1)^ = 0.01 on a possible evolu¬ 
tion of the mass variance from 0 < z < 1. Note that this 
result is robust to the assumed true value of the scatter 
when quoted as a constraint on the mass variance and 
not the rms scatter. 

However, self-calibration is not a replacement for sim¬ 
ulated catalogues, cross-calibration techniques from so- 
called direct mass measurements M, and monitoring 
scatter in observable-observable scaling relations. It is 
instead an internal consistency check on their assump¬ 


tions and the simplifying assumptions in this study. We 
have assumed that the observable-mass distribution is 
a Gaussian in InM and that its parameters depend in a 
known way on mass at a given redshift for at least a range 
in masses that is greater than the expected scatter. Fur¬ 
thermore, for low mass clusters detected optically (e.g. 
[r^ l or through leasing (e.g. 0|), the assumption of a 
one-to-one mapping of objects identihed by mass to those 
identihed by the observable breaks down since confusion 
and projection will cause many small mass objects in a 
given redshift range to be associated with a single object 
in the observable (see e.g. IlillS)- 

Without such simplifying assumptions, true self¬ 
calibration is impossible (see e.g. [3). Still, the ideas 
underlying self-calibration will be useful in revealing vio¬ 
lations of the assumed form of the distribution of cluster 
observables given the cluster masses and prevent misin¬ 
terpretation of the data. 
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